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We present a theoretical analysis of the acoustic radiation force on a single small particle, either 
a thermoviscous fluid droplet or a thermoelastic solid particle, suspended in a viscous and heat- 
conducting fluid medium. Our analysis places no restrictions on the length scales of the viscous 
and thermal boundary layer thicknesses 5^ and relative to the particle radius a, but it assumes 
the particle to be small in comparison to the acoustic wavelength A. This is the limit relevant 
to scattering of sound and ultrasound waves from micrometer-sized particles. For particles of size 
comparable to or smaller than the boundary layers, the thermoviscous theory leads to profound 
consequences for the acoustic radiation force. Not only do we predict forces orders of magnitude 
larger than expected from ideal-fluid theory, but for certain relevant choices of materials, we also 
find a sign change in the acoustic radiation force on different-sized but otherwise identical particles. 
This phenomenon may possibly be exploited in handling of submicrometer-sized particles such as 
bacteria and vira in lab-on-a-chip systems. 


I. INTRODUCTION 

The acoustic radiation force is the time-averaged force 
exerted on a particle in an acoustical field due to scatter¬ 
ing of the acoustic waves from the particle. Theoretical 
studies of the acoustic radiation force date back to King 
in 1934 [1] and Yosioka and Kawasima in 1955 [2], who 
considered the force on an incompressible and a com¬ 
pressible particle, respectively, in an inviscid ideal fluid. 
Their work was summarized and generalized in 1962 by 
Gorkov [3], however, with the analysis still limited to 
ideal fluids and valid only for particles with a radius a 
much smaller than the acoustic wavelength A. 

In subsequent work, Doinikov developed general the¬ 
oretical schemes for calculating acoustic radiation forces 
including viscous and thermoviscous effects [4-6]. The 
direct applicability of these studies is hampered by the 
generality of the developed formalism, and analytical ex¬ 
pressions are given only in the special limits of ^ ^ a ^ A 
and a S X, where S is the boundary layer thick¬ 
nesses. Similarly, the work of Danilov and Mironov, in¬ 
cluding viscous effects, only provides analytical expres¬ 
sions in these two limits [7]. However, micrometer-sized 
particles at kHz or MHz frequency relevant to acoustic 
levitation [8-11] and lab-on-a-chip applications [12-24] 
are outside these limits, because then 6 a X. This 
more general case was subsequently studied analytically 
by Settnes and Bruns including viscous boundary layers 
of arbitrary size [25]. 

In this work we extend the radiation force theory for 
droplets and elastic particles to include the effect of both 
viscosity and heat conduction, thus accounting for the 
viscous and thermal boundary layers of thickness 6^ and 
respectively, and we give closed-form analytical ex- 
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pressions in the limit of ^ with no further re¬ 

strictions between ^g, 6^, and a. Our approach to the full 
thermoviscous scattering problem follows that of Epstein 
and Carhart from 1953 [26]. The scope of their work was 
a theory for the absorption of sound in emulsions such 
as water fog in air. In 1972, Allegra and Hawley fur¬ 
ther developed the theory to include elastic solid parti¬ 
cles suspended in a fluid in order to calculate attenuation 
of sound in suspensions and emulsions [27]. The semi¬ 
nal work of these authors have become known as ECAH 
theory within the field of ultrasound characterization of 
emulsions and suspensions, and combined with the mul¬ 
tiple wave scattering theories of Refs. [28, 29] it has been 
applied to calculate homogenized complex wavenumbers 
of suspensions and emulsions [30, 31]. 

The field of ultrasound characterization driven by en¬ 
gineering applications and the field of acoustic radiation 
forces have developed in parallel with little overlap. In¬ 
deed, the scopes of the work in the two fields are very 
different. In the works of Epstein and Carhart and Alle¬ 
gra and Hawley, there is no mention of acoustic radiation 
forces [26, 27]. However, the underlying scattering prob¬ 
lem of a particle suspended in a fluid remains the same, 
and having once solved for the amplitude of the prop¬ 
agating scattered wave, the acoustic radiation force on 
the particle may be obtained from a far-held calculation. 
In the far held, the propagating scattered held changes, 
when taking into account the thermoviscous scattering 
mechanisms, including boundary layer losses and excita¬ 
tion of acoustic streaming in the vicinity of the particle. 
In this work we will elucidate this approach, as it leads 
to a particularly simple and valuable formulation for the 
acoustic radiation force in the long-wavelength limit [25]. 

Considering the success of the ECAH method to de¬ 
scribe attenuation of sound in emulsions and suspen¬ 
sions, we can with great conhdence apply the method 
to analyze the consequences of thermoviscous scattering 
on the acoustic radiation force. Nevertheless, we hnd a 
need to re-examine the problem of thermoviscous scat- 
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FIG. 1. (Color online) Sketches of the physical mechanisms responsible for various multipole components in the scattering of an 
incident acoustic wave on a particle, (a) Compressibility contrast: the incident periodic pressure field compresses the particle 
relative to the fluid, which leads to monopole radiation, (b) Thermal contrast: the incident periodic temperature field leads 
to thermal expansion of the particle relative to the fluid, which also gives rise to monopole radiation and the development of 
a diffusive thermal boundary layer (pink), (c) Density contrast: a difference in inertia between particle and fluid causes the 
particle to oscillate relative to the fluid, which gives rise to dipole radiation and the development of a viscous boundary layer 
(blue), (d) Particle resonances: acoustic wavelengths comparable to the particle size leads to complex shape changes, which 
gives rise to multipole radiation and a complex thermoviscous boundary layer (pink/blue). 


tering in order to apply the theory to the problem of 
acoustic radiation forces in a clear and consistent man¬ 
ner. One point of clarification relates to an ambiguity in 
the thermoelastic solid theory presented by Allegra and 
Hawley [27], where no clear distinction is made between 
isothermal and adiabatic solid parameters, thus tacitly 
implying 7 = c^/cy = 1 in solids. Here, we will provide 
a self-consistent treatment of thermoviscous scattering 
that clarifies this issue and allows for ease of comparison 
with existing acoustic radiation force theories. 

Before proceeding with the mathematical treatment, 
we refer the reader to Fig. 1, which illustrates the physi¬ 
cal mechanisms responsible for the monopole, dipole, and 
multipole scattering from a particle subject to a periodic 
acoustic field [31]. The final results for the acoustic radia¬ 
tion force are presented in terms of corrected expressions 
for the monopole and dipole scattering coefficients /q and 
fi. This approach allows for an easy comparison to the 
ideal fluid theory and moreover, as shown by Settnes and 
Bruns [25], it provides a simple way of evaluating acous¬ 
tic radiation forces for any given incident acoustic field. 


TABLE I. References to analytical expressions derived in this 
paper for the monopole and dipole scattering coefficients /q 
and fi in the long-wavelength limit a X. For any given 
incident acoustic field, the acoustic radiation force is 

calculated using Eq. (5) with these expressions for /q and fi. 


Size of particle and boundary layers 

fo 

fi 

Thermoviscous droplet: 

Arbitrary particle size 

Eq. (59) 

Eq. (68) 

Small-width boundary layers 

Eq. (60) 

Eq. (69) 

Zero-width boundary layers 

Eq. (61) 

Eq. (72) 

Point-particle limit 

Eq. (62) 

Eq. (73) 

Thermoelastic particle: 

Arbitrary particle size 

Eq. (64) 

Eq. (70) 

Small-width boundary layers 

Eq. (66) 

Eq. (71) 

Zero-width boundary layers 

Eq. (67) 

Eq. (72) 

Point-particle limit 

Eq. (65) 

Eq. (73) 


To this end. Table I provides an overview of the equations 
needed to evaluate the thermoviscous acoustic radiation 
force on small droplets or solid particles. 

II. BASIC CONSIDERATIONS ON THE 
ACOUSTIC RADIATION FORCE 

We consider a single particle or droplet suspended in an 
infinite, quiescent fluid medium with no net body force, 
but perturbed by a time-harmonic acoustic field with an¬ 
gular frequency cj. The density, velocity, and stress of the 
perturbed fluid is denoted p, t;, and cr, respectively. The 
region U(t) occupied by the particle, its surface 5U(t), 
and the outward-pointing surface vector n depend on 
time due to the acoustic field. The instantaneous acoustic 
radiation force is given by the surface integral of the fluid 
stress cr acting on the particle surface. However, since the 
short time scale corresponding to the oscillation period 
r is not resolved experimentally, we define the acoustic 
radiation force F^^^ in the conventional time-averaged 
sense [1-4, 7, 25], 



where the angled bracket denotes the time average over 
one oscillation period. Notice that this definition includes 
the acoustic streaming generated locally near the parti¬ 
cle, since the stresses leading to this streaming are con¬ 
tained in the fluid stress tensor cr. In contrast, by con¬ 
sidering an infinite domain, we are excluding effects of 
what Danilov and Mironov refer to as external streaming 
[7], which would be generated at the boundaries of any 
finite domain. For a given finite domain, the external 
streaming can be calculated [32], and the total force act¬ 
ing on a particle is the sum of the radiation force and the 
external-streaming-induced Stokes drag. This approach 
has been used in studies of particle trajectories and has 
been validated experimentally [33, 34]. 
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We consider a state, which is periodic in the acoustic 
oscillation period r, tantamount to requiring that any 
non-periodic phenomenon, such as particle drift, is negli¬ 
gible within one oscillation period. Usually, this require¬ 
ment is not very restrictive, as discussed in more detail 
in Section VII. For a time-periodic state, any field can be 
written as a Fourier series f(r,t) = , 

with uj = 27r/r, and the time-average of any total time 
derivative is zero, (^/(^:^)) = 0- 

A useful expression for is obtained by considering 
the momentum flux density cr — pvv entering the fluid 
volume between the particle surface and an arbi¬ 

trary static surface dCli enclosing the particle. The total 
momentum P of the fluid in this volume is the volume 
integral of pv, and because the net body force on the 
fluid is zero, the time-averaged rate of change is 


( 


dP 

dt 


a — pvv 




• n da 




a — pvv 




n da 


-irad 


( 2 ) 


Here, n is the surface vector pointing out of dfli (out of 
the fluid) and out of 9U(t) (into the fluid). The advection 
term pvv is zero at 5U(t), since there is no advection of 
momentum through the interface of the particle. Finally, 
using that the time average of the total time derivative 
is zero in the time-periodic system, we obtain 


j^rad 



a — pvv 


' n da). 


( 3 ) 


Thus, even before applying perturbation theory, the 
acoustic radiation force can be evaluated as the total mo¬ 
mentum flux through any static surface dCli enclosing the 
particle. To second order in the acoustic perturbation, 
using the expansions p = Po + Pi + P 25 v = 0 Vi V 2 ^ 
and (T = (Tq + + cr 2 , the radiation force (3) becomes 


^rad 



{(^2) - Po{viVi) 


n da. 


( 4 ) 


where we have used that the time-average of the time- 
harmonic, first-order fields is zero. 

In regions sufficiently far from acoustic boundary lay¬ 
ers, the acoustic wave is a weakly damped propagating 
acoustic mode, for which viscous and thermal effects are 
negligible. This insight was used in Ref. [25] to analyt¬ 
ically integrate Eq. (4) by placing dfti in the far field. 
In the long-wavelength limit, where the particle radius a 
is assumed much smaller than the wavelength A, i.e. for 
kQa 1 with pQ = 27r/A, it was shown that the acoustic 
radiation force may be evaluated directly from the inci¬ 
dent first-order acoustic field and the expressions for the 
monopole and dipole scattering coefficients /q and fi for 
the suspended particle, as 


F rad 3 

= —7ra 


2 k, 


Re [foPiryPin] - PoRe [fl V*n ' VVj, 


( 5 ) 


Here, and is the incident acoustic pressure and 
velocity fields evaluated at the particle position, the as¬ 
terisk denotes complex conjugation, and hCg and po are 
the isentropic compressibility and the mass density of 
the fluid medium, respectively. 

Equation (5) is valid for any incident time-harmonic, 
acoustic field, and consequently the problem of calcu¬ 
lating the radiation force on a small particle reduces to 
calculating the coefficients /q and fi. Closed, analytical 
expressions for these are given in the literature for small 
particles in the special cases of compressible particles in 
ideal fluids [2, 3] and compressible particles in viscous 
fluids [25]. Moreover, /q and fi can be extracted from 
Ref. [5, 6] for rigid spheres and liquid droplets in ther- 
moviscous fluids for the limiting cases of very thin and 
very thick boundary layers. The main result of this paper 
is the derivation of analytical expressions for /q and fi 
for a spherical thermoviscous droplet and a thermoelastic 
particle suspended in a thermoviscous fluid without re¬ 
strictions on the boundary layer thicknesses, see Table I. 
Moreover, we provide an analysis of how is affected 
by thermoviscous effects in these cases. 

Einally, we note that since /q and fi depend only on 
frequency and material parameters, expression (5) for the 
radiation force remains valid for any incident wave com¬ 
posed of plane waves at the same frequency. In the case 
of a superposition of acoustic fields 

Pj{ojj) (and similarly for i;^^) at different frequencies 
the resulting radiation force is obtained by summing over 
the forces obtained from Eq. (5) for each frequency. 


N 


F™'* = -ira y] [-^ Re [/q (wj) P*j'^Pj\ 


i=i 


-poRe [/r(wj) V* ■ Vvj 


(6) 


This generalization of Eq. (5) provides a way to evaluate 
the acoustic radiation force on a single particle regardless 
of the complexity of the incident field. 


III. THERMOVISCOUS PERTURBATION 
THEORY OF ACOUSTICS IN FLUIDS 

The starting point of the theory is the first law of ther¬ 
modynamics and the conservation of mass, momentum, 
and energy. Introducing the thermodynamic variables 
temperature T, pressure p, density p, internal energy 5 
per mass unit, entropy s per mass unit, and volume per 
mass unit 1/p, the first law of thermodynamics with s 
and p as independent variables becomes 

ds = T ds — pd[- J = Tds + dp. (7) 
\PJ p 

For acoustic wave propagation it is often convenient to 
use T and p as independent thermodynamic variables. 
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This is obtained by a Legendre transformation of the in¬ 
ternal energy 5 per unit mass to the Gibbs free energy g 
per unit mass, g = e — Ts ^ pK 

Besides the first law of thermodynamics, the governing 
equations of thermoviscous acoustics requires the intro¬ 
duction of the velocity field v and the stress tensor cr of 
the fluid. The latter can be expressed in terms of v, p, 
the dynamic shear viscosity the bulk viscosity 7 ^^, and 
the viscosity ratio /3 = p^/p-hl/3, as 


cr = —p I r, (8a) 

T = T] Vv + (Vi;)^ + (/5 — 1)^ (V • v) I. (8b) 


s = s{p^ T). The total differentials of pi and Si are 

which may be linearized so that the partial derivatives of 
p and s refer to the unperturbed state of the fluid. This 
leads to the introduction of the isothermal compressibil¬ 
ity the isobaric thermal expansion coefficient and 
the specific heat capacity at constant pressure c^. 


Here, I is the unit tensor and the superscript ”T” indi¬ 
cates tensor transposition. The tensor r is the viscous 
part of the stress tensor assuming a Newtonian fluid [35]. 

Considering the fluxes of mass, momentum and energy 
into a small test volume, we use Gauss’s theorem to for¬ 
mulate the general governing equations for conservation 
of mass, momentum and energy in the fluid under the 
assumption of no net body forces and no heat sources, 

dtP = '^ ■ [-pv], (9a) 

dtipv) = V • [<7 - pvv], (9b) 

dt{ps + ^pv"^) = V • \y ■ a + - p{e + \v^)v\. 

(9c) 


ds \ 

(13) 

Moreover, {ds/dp)j- = —a^jp^ which may be derived as 
a Maxwell relation differentiating g after p and T. Thus, 
the linearized form of Eq. (12) is 

Pi = po i^T Pi ~ Po Ti, (14a) 

Cr, .X 

Si = t^T^-^Pi^ (14b) 

To Po 

We further introduce the isentropic compressibility 
and the specific heat capacity at constant volume Cy, 


Kj'P' — 


I f dp 


P \dp 


ar. 


P\dT), 


c =T 

1 '-'p 


Here, we have introduced the thermal conductivity 
assuming the usual linear form for the heat flux given by 
Fourier’s law of heat conduction. 


A. First-order equations for fluids 

The zeroth-order state of the fluid is quiescent, homo¬ 
geneous, and isotropic. Then, treating the acoustic field 
as a perturbation of this state in the acoustic perturba¬ 
tion parameter given by 

eac = — « 1, (10) 

Po 

we expand all fields as p = po + with Vq = 0. 

The zeroth-order terms drop out of the governing equa¬ 
tions, while the first-order mass, momentum, and energy 
equations obtained from Eqs. (7) and (9) become 

dtPi =-Po^-Vi: (11a) 

PodtVi = -Vpi + (V • t;/ , (11b) 

Po^odfSi = (11c) 

It will prove useful to eliminate the variables p^, p^, and 
to end up with only two equations for the variables 
Vi and Ti. To this end, we combine Eq. (11) with the 
two thermodynamic equations of state p = p(p, T) and 


-Kl)/ 



(15a) 


Then the following two well-known thermodynamic iden¬ 
tities may be derived [36], 




7 = —= 1 + 
Cy 


2 rri 
dpi Q 

PoCpt^s 


(16) 


To proceed with the reduction of Eq. (11), we first dif¬ 
ferentiate Eq. (11b) with respect to time and substitute 
= V(V • Vi) — V X V X Then Eq. (14) is used 
to eliminate pi and Si in Eqs. (11b) and (lie), followed 
by elimination of d^pi using Eq. (11a). The resulting 
equations for Vi and are 


dt 


+ (l + /3)z/oa, V(V-i;i) 


PO^T 

d 

+ i^odf^ X V X Vi = -(17a) 

Po^T 


7AhV"Ti - dy, = 

dr. 


(17b) 


where we have introduced the momentum diffusion con¬ 
stant Vo ^^4 the thermal diffusion constant 


^0 — —: ^th — -• 

Po PoCp 


( 18 ) 
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B. Potential equations for fluids 


The velocity field is decomposed into the gradient 
of a scalar potential (j) (the longitudinal component) and 
the rotation of a divergence-free vector potential -0 (the 
transverse component), 

Vi = + V X -0, with V • -0 = 0. (19) 


Inserting this well-known Helmholtz decomposition into 
Eq. (17a) leads to the equation 


dU- 


1 


= V X 


-h (1 + ) V 0 H- —dfTi 

PohvT 


-dt'ip + vodtV'^ip 


( 20 ) 


In general, both sides of the equation must vanish sep¬ 
arately, which leads to two equations. Combining these 
with Eq. (17b), into which Eq. (19) is inserted, leads to 
the following form of Eq. (17), 

dU = (+ (1 + l3Hdt] VV - (21a) 

VPO^T / PO^T 

(21b) 

(Xp 

dt'ip = i/oVV- (21c) 

In the adiabatic limit, for which = 0, the well-known 
adiabatic wave equation for 0 is obtained by inserting 
Eq. (21b) into (21a), from which the adiabatic speed of 
sound c for longitudinal waves is deduced, 

r; = / . (22) 

V Po^s 

In the isothermal case, for which =0, the wave 
equation (21a) instead describes waves traveling at the 
isothermal speed of sound c/^Py = Eor ultra¬ 

sound acoustics, sound propagation in the bulk of a fluid 
is generally very close to being adiabatic. 


contrast to the per-mass quantities given by lowercase 
letters. The first law of thermodynamics reads 


dS = TdS (T^jdu^j^ (23) 

where £ is the internal energy per unit volume, S is the 
entropy per unit volume, and T is the temperature. The 
work done by the internal stresses per unit volume is 
equal to —where we have introduced the strain 
tensor which for small displacements is given by 

Uij = ^ [diUj + djUi] . (24) 

Transforming the internal energy per unit volume £ to 
the Helmholtz free energy per unit volume F = £ — TS, 
where temperature T and strain u^j are the independent 
variables, the first law becomes dF = —SdT + a^jdu^j. 

Consider the undeformed state of an isotropic, ther¬ 
moelastic solid at temperature Tq in the absence of ex¬ 
ternal forces. The free energy F is then given as an ex¬ 
pansion in powers of the temperature difference T — Tq 
and the strain tensor u^j. To linear order, the stress ten¬ 
sor a^j = and the entropy S = — (f^)^ become 


ij T ^ c^^^kk^ij 5 

(25a) 

S{T) = SoiT) + (25b) 


ap{T — Tq) ^ 


E 


where Sq{T) is the entropy of the undeformed state at 
temperature T, while E and a are the isothermal Young’s 
modulus and Poisson’s ratio, respectively. The isother¬ 
mal compressibility Kj- of the solid is given in terms of E 
and (7 as. 


_ 3(1 -2cr) 
“ E 


(26) 


A. Linear equations for solids 


IV. THERMOELASTIC THEORY OF 
ACOUSTICS IN ISOTROPIC SOLIDS 

A thermoelastic solid may be deformed by the action 
of applied forces or on account of thermal expansion. Fol¬ 
lowing Landau and Lifshitz [37], we describe the defor¬ 
mation of a solid elastic body using the displacement field 
n, which describes the displacement u{r^t) of a solid el¬ 
ement away from its initial, undeformed position r to 
its new temporary position r Any displace¬ 

ment away from equilibrium gives rise to internal stresses 
tending to return the body to equilibrium. These forces 
are described using the stress tensor cr, which leads to 
the force density V • cr. In the description of the ther¬ 
modynamics of solids, it is advantageous to work with 
per-volume quantities denoted by uppercase letters, in 


In elastic solids advection of momentum and heat can¬ 
not occur, so the momentum equation in the absence of 
body forces takes the linear form pd^u = V • cr. As¬ 
suming the material parameters F, and a to be 

constant, it becomes 


E 

pd^u = -^VT + —^ 
2(1 F cr) 

a 




1 -2cr 


V(V-m) 


^VT+CtV w+(cl-Ct)V(V--u), (27) 

pKx 


where we have introduced the isothermal speed of sound 
of longitudinal waves Cl and of transverse waves Cx, 


2 _ (1 - 0-) E 2 _ 1 E 

(1 + C7)(l - 2 c 7) p ’ 2{l + a)p' 


(28a) 
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Using the decomposition n = + '^l the trans¬ 
verse and longitudinal displacements Ul with 

V • Ux = 0 and V x Ul = 0? respectively, it immedi¬ 
ately follows from Eq. (27) that in the isothermal case, 
transverse and longitudinal waves travel at the speed Cx 
and Cx, respectively. Combining Eqs. (26) and (28a) one 
obtains an important relation connecting the isothermal 
compressibility of the solid to the isothermal sound 
speeds Cx and Cx, 


1 



(28b) 


Turning to the energy equation, the amount of heat ab¬ 
sorbed per unit time per unit volume is T{dfS). If there 
are no heat sources in the bulk, the rate of heat absorbed 
is given by the influx of heat by conduction, and 

the heat equation thus becomes 


T{dtS) = -V • [-fcthVT] = (29) 

where the heat conductivity is taken to be constant. 
We rewrite this equation using expression (25b) for the 
entropy, and using that the time-derivative of Sq may be 
written as 


dt \dT )y dt T dt' ^ ’ 

where the heat capacity Cy per unit volume at constant 
volume enters through the relation Cy = T {dSQ/dT)y 
with the derivative taken for the undeformed state at 
constant volume, that is for Uf.f. = V • n = 0. Combining 
these considerations with the identity for 7 equivalent to 
Eq. (16), the heat equation (29) becomes 

Cydy + ■ u = ( 31 ) 

dp 

Einally, having eliminated all extensive thermodynamic 
variables, we return to per-mass quantities, such as Cy = 
Cy /p, and thus arrive at the coupled equations for ther¬ 
moelastic solids, 

d^Ui — CxV(V • Ui) + CxV X V X Ui = - —VTi, 

PqK,t 

(32a) 

lAhV^Ti - dy, = ■ Uy (32b) 

dp 

with 7 and defined in Eqs. (16) and (18), and the 
linearity emphasized by the addition of subscripts ” 1 ” 
to the field variables. In this form, the thermoelastic 
equations (32) correspond to the fluid equations (17). 


fluid case, we make a Helmholtz decomposition of this 
velocity field in terms of the velocity potentials (j) and -0 

dfUi = V0 + V X -0, with V • -0 = 0. (33) 

Inserting this into Eq. (32) and following the procedure 
leading to Eq. (21) for fluids, we obtain the corresponding 
three equations for solids. 


ptlry 

(34a) 

dy, = TAhV'Ti - 

(34b) 

dp 


dlxl, = 

(34c) 


The main difference between the fluid and the solid case 
is in Eq. (34c) for the vector potential -0, which now takes 
the form of a wave equation describing transverse waves 
traveling at the transverse speed of sound Cx instead of 
the diffusion equation ( 21 c). 

The usual adiabatic wave equation for the scalar po¬ 
tential (j) is obtained in the limit of = 0 combining 
Eqs. (34a) and (34b), and the speed c of adiabatic, lon¬ 
gitudinal wave propagation in an elastic solid becomes 


2 2 , 

c — Cx + 


7-1 

PO^T 


(35) 


Eor most solids 7 — 1 <C 1, leading to a negligible differ¬ 
ence between the isothermal Cx and the adiabatic c, the 
latter being closest to the actual speed of sound measured 
in ultrasonic experiments. 


V. UNIFIED POTENTIAL THEORY OF 
ACOUSTICS IN FLUIDS AND SOLIDS 


The similarity between the potential equations (21) 
and (34), allows us to write down a unified potential the¬ 
ory of acoustics in thermoviscous fluids and thermoelastic 
solids. The main result of this section is the derivation 
of three wave equations with three distinct wavenum¬ 
bers corresponding to three modes of wave propagation, 
namely two longitudinal modes describing propagating 
compressional waves and damped thermal waves, respec¬ 
tively, and one transverse mode describing a shear wave, 
which is damped in a fluid but propagating in a solid. 

We work with the first-order fields in the frequency 
domain considering a single frequency cj. Using complex 
notation, we write any first-order field gi{r,t) as 


9 i{r,t) = gj{r)e (36) 

B. Potential equations for solids 

Assuming this form of time-harmonic first-order fields. 
The time derivative d^Ui of the displacement field Ui Eqs. (21a) and (34a) lead to expressions for the temper- 

describes the velocity field in the solid. Analogous to the ature field in a fluid (fl) and a solid (si), respectively. 
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in terms of the corresponding scalar potential 


fl _ lUJpol^T 

- 1-1 — - 

Clrn 


si _ 

-1 ~ —I- 


c i-i7r 2 

—-V (j 

u 7 


.+4v^ 

U! 


(37a) 

(37b) 


Here, we have introduced the dimensionless bulk damp¬ 
ing factor Fg accounting for viscous dissipation in the 
fluid. For convenience, we also introduce the thermal 
damping factor F^ accounting for dissipation due to heat 
conduction both in fluids and in solids. These two bulk 
damping factors are given by 


F 


s 


(1 + l3)l'oUJ 



(38) 


Substituting expression (37a) for Tf into Eq. (21b), 
or expression (37b) for Ti into Eq. (34b), and assuming 
time-harmonic fields Eq. (36), we eliminate the tempera¬ 
ture field and obtain a bi-harmonic equation for the scalar 
potential 0 , 


= O 5 with = —^ (39a) 
c 

where we have introduced the undamped adiabatic 
wavenumber /cq = and where the parameters o^xi 

and /^xi for fluids (xl = fl) and solids (xl = si) are 


Qffl — -i(l - i7Fs)Ft, — 1 - i(Fs -|- qF^), (39b) 

^si = -i(l + X)Ft, = 1 - iTFf (39c) 


Here, we have used the relation (35) for solids and further 
introduced the parameters X and y. 


X=(7-l)(l-x) 

1 


X = 


= 1 - 


PqKjqC 


4 Cx 


(39d) 

(39e) 


the latter equality following from combining Eq. (35) 
with Eq. (28b) and using = 7/^5 from Eq. (16). Note 
that for fluids y = 1 , = 0 , and X = 0 . 

The bi-harmonic equation (39a) is factorized and writ¬ 
ten on the equivalent form 


(V^ + kl){V^ + kl)<j> = 0, (40a) 


and thus the wavenumbers and k^ are obtained from 

resulting in 


kc — ^xi + (/^xi ~ 4(axi)^^^ 

kl = 2kl ki - {dl, - 


(40b) 

(40c) 


with ’’xl” being either ”fl” for fluids or ”sl” for solids. 

In the frequency domain, the equation for the vector 
potential -0, Eq. (21c) for fluids and Eq. (34c) for solids. 


2 2 

can be written as V ^l) ^ k^ if) = 0, which describes a 
transverse shear mode with shear wavenumber k^. By 
introducing a shear constant 770 ? which for a fluid is the 
dynamic viscosity, and for a solid is defined as, 

2 

770 = i (solid), (41a) 

UJ 

the shear wavenumber k^ is given by the same expression 
for both fluids and solids, 

kl = —— (fluid and solid). (41b) 

Vo 


A. Wave equations and modes 

The general solution cj) of the bi-harmonic equa¬ 
tion (40a) is the sum 


4 = 4c + 4t (42) 

of the two potentials (j)^ and 0^, which satisfy the har¬ 
monic equations 


+ kcdc = 0, (43a) 

V7t + kUt = 0, (43b) 

where describes a compressional propagating mode 
with wavenumber k^, while describes a thermal mode 
with wavenumber k^. These two scalar wave equations 
together with the vector wave equation for describing 
the shear mode with wavenumber /Cg, 

+ /Cg-0 = 0 , (43c) 

comprise the full set of first-order equations in poten¬ 

tial theory. These wave equations, coupled through the 
boundary conditions, govern acoustics in thermoviscous 
fluids and thermoelastic solids. The distinction between 
fluids and solids is to be found solely in the wavenumbers 
of the three modes. 


1 . Approximate wavenumbers for fluids 


Eor most systems of interest, Fs,F^ ^ 1 allowing a 
simplification of the expressions for k^ and k^ in Eq. (40). 
To first order in Fg and F^ one finds 


kc 

kt 


UJ 

c 


i + ^[rs + (7-i)rt] 


(l+i) 


i + ^( 7 -i)(r,- 


(1+i) 



(44a) 

(44b) 

(44c) 
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where we have introduced the thermal diffusion length 
and the momentum diffusion length S^. Heat and mo¬ 
mentum diffuses from boundaries, such that the charac¬ 
teristic thicknesses of the thermal and viscous boundary 
layers are and respectively, given by 



For water at room temperature and 2 MHz frequency, 
^ 0.4 pm, (5^ 0.2 pm, and A 760 pm. Con¬ 

sequently, the length scales of the thermal and viscous 
boundary layer thicknesses are the same order of mag¬ 
nitude and much smaller than the acoustic wavelength. 


With ko = 

uj/ewe 

note that 




rs = 

2 + /^)(^0^s)^5 

Tt 

= ^{koSif, 

(46) 

and consequently 






Ts' 

(ko^af ~ 

K 

2 

«1, 

(47a) 




kf. 

2 

« 1. 

(47b) 


In the long-wavelength limit of the scattering theory to 
be developed, we expand to first order in /cq^s 
and thus neglect the second-order quantities Fg and F^. 
For water at room temperature and MHz frequency one 
finds /cq^s ^ ^ 10 ~^, and Fg ^ F^ ^ 10~^. 

Clearly, the compressional mode with wavenumber 
kc describes a weakly damped propagating wave with 
Im[/Cc] ^ ujc. In contrast, Im[/ct] — Re[/ct] 

for the thermal mode and Im[/Cs] = Re[/Cs] for the shear 
mode, which correspond to waves that are damped within 
their respective wavelengths. Hence, these modes de¬ 
scribe boundary layers near interfaces of walls and par¬ 
ticles, which decay exponentially away from these inter¬ 
faces on the length scales set by 6^ and 6^. 


B. Acoustic fields from potentials 


For a given thermoacoustic problem, the boundary 
conditions are imposed on the acoustic fields Vi, , and 
a I and not directly on the potentials We 

therefore need expressions for the acoustic fields in terms 
of the potentials in order to derive the boundary condi¬ 
tions for the latter. 

The velocity fields follow trivially from the Helmholtz 
decompositions and are obtained from the same expres¬ 
sion in both fiuids and solids 


I’l = V+ 4) + V X V’, (49) 

where Vi = —iujUi for solids. 

A single expression for in terms of (j)^ and (j )^, valid for 
both fiuids and solids, is obtained from Eq. (37) in combi¬ 
nation with Eqs. (40) - (43) by introducing the material- 
dependent parameters and 6^, 


T 

_ i^(7 - 1) , _ 1 

— 2 ’ • 

apC X<^pAh 


(50a) 

(50b) 


Here, we have neglected Fg and F^ relative to unity. Note 
that the ratio b^/b^ ^ 1. 

In a fiuid, the pressure field pi is obtained by inserting 
Eq. (19) into the momentum equation (11b) and using 
the wave equations (43), 


Pi = iw/?o(0c + 4>t) - (1 + l3)vo{kc4>c + (51) 


Inserting this expression into Eq. (8a), the stress tensor 
for fiuids becomes. 


<7i = 7?o (2fc^ - + (2fct - fcg )0t 

+ Vo [Vvi + (Vvi)"^ 


(52) 


2. Approximate wavenumbers for solids 


Similar to the fiuid case, we use the smallness of the 
thermal damping factor, F^ ^ 1, to expand the exact 
wavenumbers of Eq. (40). To first order we obtain 


UJ 

c 


1 + 2(9'- i)xrt 


k,= 


(1 + i) 


<5t 


1 7 Ft 

1 + -- * 


8(1-X) 


fc, = -. 

Ct 


(48a) 

(48b) 

(48c) 


An important distinction between a fiuid and a solid is 
that a solid allows for propagating transverse waves while 
a fiuid does not. This is evident from the shear mode 
wavenumber /Cg, which for solids is purely real, 
while for fiuids Im[/Cs] = Re[/Cs] = l/c^g. 


where can be expressed by the potentials through 
Eq. (49). This expression also holds true for the solid 
stress tensor Eq. (25a) using the shear constant tjq, 
Eq. (41a), and the velocity field Vi = —iujUi, Eq. (33). 
This conclusion is obtained by inserting Eq. (37b) for T® 
into Eq. (25a) for cr^ and using the wave equations (43). 


VI. SCATTERING FROM A SPHERE 

The potential theory allows us in a unified manner to 
treat linear scattering of an acoustic wave on a spherical 
particle, consisting of either a thermoelastic solid or a 
thermoviscous fiuid. The system of equations describing 
the general case of an arbitrary particle size is given, and 
analytical solutions are provided in the long-wavelength 
limit a,(5g,(5t ^ A. In this limit, the particle and bound¬ 
ary layers are much smaller than the acoustic wavelength, 
but the ratios d^/a and 6^/a are unrestricted. This is 
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Fluid medium / 


Propagating waves 

(compressional mode) 
Incident wave: 
Reflected wave: 



FIG. 2 . (Color online) A compressional wave propagating in a thermoviscous fluid medium with parameters 
Po?'^ 0 ? CKp? Cp, 7 , and /cth? is incident on a thermoviscous fluid droplet with parameters po, ? 7 o, Op, Cp, 7 ^, and /cthj which 
results in a compressional scattered wave and highly damped thermal and shear waves both outside in the fluid medium 
(0r, "^s) cind inside in the fluid droplet (0c, Viscous and thermal boundary layers are described by the highly damped 

waves both outside and inside the fluid droplet. In the long-wavelength limit the droplet radius a and the boundary layer 
thicknesses ^ire mutually unrestricted, but all much smaller than the acoustic wavelength A. For a thermoelastic 

particle, the shear mode 'ipg describes a propagating transverse wave instead of an internal viscous boundary layer. 


essential for applying our results to micro- and nanopar¬ 
ticle acoustophoresis. In particular, we derive analytical 
expressions for the monopole and dipole scattering coeffi¬ 
cients /o and fi , which together with the incident acous¬ 
tic field serve to calculate the acoustic radiation force as 
shown in Section II and summarized in Table I. 


A. System setup 

We place the spherical particle of radius a at the center 
of the coordinate system and use spherical coordinates 
(r, 0^ ip) with the radial distance r, the polar angle 0^ and 
the azimuthal angle (p. We let unprimed variables and 
parameters characterize the region of the fluid medium, 
r > while primed variables and parameters character¬ 
ize the region of the particle, r < a. For example, the 
parameter K,'g is the compressibility of the particle, while 
Kg is the compressibility of the fluid medium. Ratios of 
particle and fluid parameters are denoted by a tilde, e.g. 

Due to linearity, we can without loss of 
generality assume that in the vicinity of the particle, the 
incident wave is a plane wave propagating in the positive 
^-direction, 0^ = 0o = 0o The fields do 

not depend on cp due to azimuthal symmetry. 

B. Partial wave expansion 

The solution to the scalar and the vector wave equa¬ 
tions Eq. (43) with wavenumbers /c, Eqs. (44) and (48), 
in spherical coordinates is standard textbook material. 
Avoiding singular solutions at r = 0 and considering 
outgoing scattered waves, the solution is written in terms 
of spherical Bessel functions j^{kr), outgoing spherical 
Hankel functions hj^{kr)^ and Legendre polynomials 
P^(cos6>). As a consequence of azimuthal symmetry. 


only the (^-component of the vector potential is non-zero, 
'?0(r) = '0s(r, ^) e^. The solution is written as a partial 
wave expansion of the incident propagating wave 0^, the 
scattered reflected propagating wave 0^., the scattered 
thermal wave 0^, and the scattered shear wave 

In the fluid medium, r > a 


= + I)in(^c0-Pn{cos6»), (53a) 

n=0 

oo 

(l>v = <f>oYl + ^)AnK{Kr)Pn{cose), (53b) 

n=0 

oo 

0t = 00 E + ^)BnK{K'r)Pnicose), (53c) 

n=0 

oo 

0s = 00 E + l)C„hn{k^r)d0P„{cosd), (53d) 

n=0 

In the particle, r < a 


0c = 00 E + l)Xin(fcc?’)-Pn(cOS0), (53e) 

n=0 

oo 

0t = 00 E + ^)Bnjn{Kr)Pn{cosO), (53f) 

n=0 

oo 

0G 00 E + l)Cnin(^k)0e-Pn(cOS0), (53g) 

n=0 

where the parameter 0o is an arbitrary amplitude of the 
incident wave with unit m^s~^. The different components 
of the resulting acoustic field are illustrated in Eig. 2. 
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C. Boundary conditions 


2 _ 2 / 

d 17^ 7^ 17^ 7^ 


Neglecting surface tension, the appropriate boundary 
conditions at the particle surface are continuity of veloc¬ 
ity, normal stress, temperature, and heat flux. Assum¬ 
ing sufficiently small oscillations, see Section VIIC, the 
boundary conditions are imposed at r = a. 


Vo 


- 2xl)j^{x^) - 2 xIj^{x^) 

+ ^nVo (^s - - 2x\}i'^{x^) 

+ BnVo [(^S - 2xt)h^(xt) - 2xth^(xt) 


V^r ^10 ^10 •) ^1 (54a) 

^Irr ^Irri ^19r ^lOr") ^th^r^l (54b) 

The boundary conditions are expressed in terms of the 
potentials using Eqs. (49), (50) and (52). The compo¬ 
nents of velocity and stress in spherical coordinates are 
given in Appendix A. 

It is convenient to introduce the non-dimensionalized 
wavenumbers and x^ for the medium, and x[, 

and Xg for the particle, 

x^ = Xi = /c^a, x^ = (55a) 

x'^ = /c^a, x[ = /c(a, Xg = k'^a. (55b) 

Inserting the expansion (53) into the boundary condi¬ 
tions (54), and making use of the Legendre equation 
(Cl), we obtain the following system of coupled linear 
equations for the expansion coefficients in each order n, 

avi^ = av'ij. 

Xcj'n{Xc) + AnXch'„{Xc) + BnXthn{x^,)-C„n{n+l)hn{Xs) 

= J^nXcfniXc) + “ C'n«(n+l)j„(a;'), (56a) 

avi0 = av'iQ 

jn{Xc') “t“ ^n^n(^c) “t“ -^n^n{^t) [^s^n(^s) ^n(^s)] 

-^nJn(Xc) [^sJn{^s) Jn(^s)] 5 

(56b) 

Ti = t; 


+ 2n(n + l)CnT}o [xs^ix^) - h^ix^)] 

= Kvo (Xs - ‘2x'^)jniXc) - ‘^x'^fnix'c) 

+ B'nVo {x? - 2a;f )j„(a;() - 2a;f j"(a;() 
+ 2n(n + l)Cnrio [x'jn{xi) - j„(a;')] • 


(56f) 


Here, primes on spherical Bessel and Hankel functions 
indicate derivatives with respect to the argument. The 
equations are valid for both a fluid and solid particle, 
with 77o being the viscosity for a fluid particle and the 
shear constant Eq. (41a) for a solid particle. 

Eor n = 0, the boundary conditions for Viq and aiQ^ are 
trivially satisfied because there is no angular dependence 
in the zeroth-order Legendre polynomial, Po(cos^) = 1. 
Consequently, ^/^g = 0, and we are left with four equations 
with four unknowns, namely Eqs. (56a), (56c), (56d), and 
(56f) with Cq = Cq = 0. 

The linear system of equations (56) may be solved for 
each order n yielding the scattered field with increas¬ 
ing accuracy as higher-order multipoles are taken into 
account, an approach referred to within the field of ul¬ 
trasound characterization of emulsions and suspensions 
as ECAH theory after Epstein and Carhart [26] and Al- 
legra and Hawley [27]. However, care must be taken due 
to the system matrix often being ill-conditioned [38]. 

The long-wavelength limit is characterized by the small 
dimensionless parameter 5, given by 

5 = k^a = 27r Y < 1, (57) 

A 


^cJn(^c) ^n^c^n(^c) -^n^t^n(^t) 

= KKjniXc) + Bnb[j„{x[), (56c) 

gfcth^r = aKhdrTj 


KhKXcj'niXc) + A^KhbcXch’niXc) + 

= KKhbcx'cjhiXc) + Bnk[i,b[x[jn{x0, (56d) 

2 2 / 
a (JlQrp (2 (Jigrp 


Vo [Xcfni^c) - jn{Xc)] + A^VO “ ^n(^c)] 

4" ^nVo [^t^n(^t) ^n(^t)] 


- ^C^Vo [xsKM ^ {n -hn-2)h^(xg) 

AfiVo [^c4n(^c) 4n(^c)] T ^nVo 4n(^t)] 


- 2 in(^s) ^n- 2)j^(xg)J , (56e) 


In this limit, the dominant contributions to the scattered 
field are due to the n = 0 monoj)ole and the n = 1 dipole 
terms, both proportional to 5'^, while the contribution 
of the nth-order multipole for n > 1 is proportional to 


D. Monopole scattering coefficient 

To obtain the monopole scattering coefficient /q in 
Eq. (5), we solve for the expansion coefficient Aq in 
Eq. (56) and use the identity /q = 3i x^Aq. The f^- 
coefficients are traditionally used in work on the acoustic 
radiation force, while the A^-coefficients are used in gen¬ 
eral scattering theory. 

The solution to the inhomogeneous system of linear 
equations for n = 0 involves straightforward but lengthy 
algebra presented in Appendix B 1. In Eq. (B8) is given 
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the general analytical expression for /q in the long- 
wavelength limit valid for any particle. In the follow¬ 
ing, this expression is given in explicit, simplified, closed 
analytical form for a thermoviscous droplet and a ther¬ 
moelastic particle, respectively. 


1. A thermoviscous droplet in a fluid 


For a thermoviscous droplet in a fluid in the long- 
wavelength limit, the particle radius a and the viscous 
and thermal boundary layers both inside ((^g, and out¬ 
side (^s, fhe fluid droplet are all much smaller than the 
acoustic wavelength A, while nothing is assumed about 
the relative magnitudes of a, 5^^ and S[. Thus, us¬ 

ing the non-dimensionalized wavenumbers Eq. (55) and 
e = the long-wavelength limit is defined as 


1 12 1 / 1 2 2^.-1 

\Xc\ , \xc\ -- e <C 1 

and 

(58a) 

1 12 1 / 12 2^.1 
kcl , kcl ^ £ « ks 


(58b) 

which implies 



p p l^cl 

S, t, 

iTTi ^ ^ 1- 

\h\ 

(58c) 


To first order in 5 , the analytical result for the monopole 
scattering coefficient /q obtained from Eq. (B 8 ) is most 
conveniently written as 


/« = 1 - K, + 3(7 - 1) 1 - 


H{x^ , (Tt) 2 


Po Cp 

1 tan x- 


1 — ixt fan — x^_ 


H{x^,x[), (59a) 
(59b) 


-1 


where H{x^^x[) is a function of the particle radius a 
through the non-dimensionalized thermal wavenumbers 
x^ and x[. Epstein and Carhart obtained a correspond¬ 
ing result for Aq but with a sign-error in the thermal cor¬ 
rection term [26], while the result of Allegra and Hawley 
[27] is in agreement with what we present here. The fac¬ 
tor (7 — 1 ) quantifies the coupling between heat and the 
mechanical pressure waves. This factor is multiplied by 
[l—dp/(poCp)] , where the quantity = a^/{poCp), with 
unit m^/J, may be interpreted as an isobaric expansion 
coefficient per added heat unit. The thermal correction 
can only be non-zero if there is a contrast ^p 7 ^ 1 in this 
parameter. 

In the weak dissipative limit of small boundary layers 
the function H {x ^, ) is expanded to first order in S^/a 
and and using tan(x[) i, we obtain 


fS 


3 (l+i)(7-l) 
2 1 + Dipt 


cy 


th '^th ^ 


a ' 


(60a) 




(Small-width boundary layers), 

(60b) 


In the limit of zero boundary-layer thickness 5^/a 0, 

the thermal correction vanishes, and we obtain 

/o = I — kg, (Zero-width boundary layers), ( 6 Ia) 

kcP « 1 and |Xs|^, 1x7, ktl^, ^7 oo, (61b) 

which is the well-known result for a compressible sphere 
in an ideal [3] or a viscous [25] fluid. 

In the opposite limit of a point particle, a/6^^ a/6[ 0, 

we find = -{l/3)poCp, yielding 

fo = 'i- - Ks - il - l)PoCp (1 - , (62a) 

(Point-particle limit), 

|Xc|t |x 7 < 1 and |x 7 , ktit 141^ 0 . (62b) 

Since 7 > I, the correction from thermal effects in the 
point-particle limit is negative. This implies that the 
thermal correction enhances the magnitude of /q for 
acoustically soft particles (kg > I), while it diminishes 
the magnitude and eventually may reverse the sign of /q 
for acoustically hard particles (kg < I). 

Importantly, an inspection of the point-particle limit 
Eq. (62) leads to two noteworthy conclusions not pre¬ 
viously discussed in the literature. Eirstly, the thermal 
contribution to /q allows for a sign change of the acoustic 
radiation force for different-sized but otherwise identical 
particles. Secondly, the thermal contribution may re¬ 
sult in forces orders of magnitude larger than expected 
from both ideal [3] and viscous [25] theory. Eor example, 
Po ^ 1 for particles or droplets in gases leads to a ther¬ 
mal contribution to /q two orders of magnitude larger 
than I — kg. These predictions are discussed in more 
detail in Section VIII. 


2. A thermoelastic particle in a fluid 


Eor a thermoelastic particle in a fluid, the long- 
wavelength limit differs from that of a thermoviscous 
droplet Eq. (58) by the shear mode describing a prop¬ 
agating wave and not a viscous boundary layer. The 
wavelength of this transverse shear wave is comparable 
to that of the longitudinal compressional wave, and in 
the long-wavelength limit both are assumed to be large. 


kcl^, k7, |a;7 ~ < 1 and 

7iMx7, 


which implies 


1 17 2 

l-sAt? |~ 15 17 15 17 I ^ I • 

l%l l^^tl l^^tl 


(63a) 

(63b) 


(63c) 


To first order in 5 , the result Eq. (B8) for may be 
simplified as outlined in Appendix B, and one obtains 
after some manipulation 
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1 — kg 3(7 — 1) 


fo 


(y.rn 


Po^p 


X o^p 

pQCp 


A I ~ ~ /2 

4 X Crp 


1 - 


cy 


PoCphZ^ 




l + iiy-l)^;^^■H{x^,X^) 
PoCp c 


(64) 


where the function is still given by the expression in Eq. (59b) with x[ being the non-dimensionalized thermal 

wavenumber in the solid particle obtained from Eq. (48b). In the limit of a point particle, a/6^,a/6[ 0, we find 


1 - ^ - 


/o' = 


(7 - l)PoCp 
1-X' 


I — 


PoCp 


1 - 


X k^p 

PoCp 


4 X ^pks Cx 

3 c„ 


1 - 


Oirn 


PoCph 


1 - 


4 7 — 1 X dp Cx 


(Point-particle limit), (65a) 


^cl^ l^cl^ < 1 and |xt|^ \x[\^ 0. (65b) 


2 I /|2 


Remarkably, in the point-particle limit and /q dif¬ 
fer in general. However, as expected, letting Cx ^ 0 in 
Eq. (64), Jq reduces to /q Eq. (59) for all particle sizes. 

In the weak dissipative limit of small boundary layers, 
(5t, ^ a, the second term in the denominator of Eq. (64) 

is small for typical material parameters. An expansion 
in d^/a and d[/a then yields in analogy with Eq. (60), 


/o' = 1 - Ks 


3 (i+i)(7-i) / _ ^ 

2i + (i-xy/2777hA pocp) a 


(66a) 

(Small-width boundary layers), 

kcit \Xc?, < 1 < ksit ktit \xt\'^, (66b) 


simplified using Eq. (39e). In the limit S^/a 0, the 
thermal correction terms vanishes, 

/q = 1 — kg^ (Zero-width boundary layers), (67a) 
kcit WX, < 1 and |a;s|^ |a;t|^ \x[\^ oo, (67b) 


In this limit, where boundary layer effects are negligible, 
and /o are identical and, as expected, equal to the 
ideal [3] and viscous [25] results. 


E. Dipole scattering coefficient 

To obtain the dipole scattering coefficient fi in Eq. (5), 
we solve for the expansion coefficient Ai in Eq. (56) and 
use the identity fi = —dix^ Ai. In the long-wavelength 
limit, the terms involving the coefficients Bi and B[ are 
neglected to first order in 5. This reduces the system of 
equations (56) for n = I from six to four equations with 
the unknowns A ^, A'l , Ci , and C [. In Appendix B 2 we 
solve explicitly for A^. Physically, the smallness of the 
Bi- and -terms means that thermal effects are negli¬ 
gible compared to viscous effects. This is consistent with 


the dipole mode describing the center-of-mass oscillations 
of the undeformed particle. 


1. A thermoviscous droplet in a fluid 


The analytical expression for Ai in the long- 
wavelength limit for a thermoviscous droplet in a fluid, as 
defined in Eq. (58), is given in Eq. (B23) of Appendix B 2. 
This expression for Ai was also obtained by Allegra and 
Hawley [27] and, with a minor misprint, by Epstein and 
Carhart [26] in their studies of sound attenuation in emul¬ 
sions and suspensions. We write the result for the dipole 
scattering coefficient fi on a form more suitable for com¬ 
parison to the theory of acoustic radiation forces as pre¬ 
sented by Gorkov [3] and Settnes and Bruus [25], 


•fl ^ 2 (po - 1) (1 + F{xs, xj) - G(x^)) 

{2pQ + 1) [l + F(a;g,a;s)] — 3G{Xs) 


G{x,) 


x<. 



(68b) 


F{xs,xi) 


2(1 — Vo) + 


1 - ix^ 

PpXg (tanxg -Xg) 

(3 — x'g) tanXg — 3Xg 


(68c) 


Even though no thermal effects are present in /f, Eq. (68) 
is nevertheless an extension of the result by Settnes and 
Bruus [25], since we have taken into account a finite vis¬ 
cosity in the droplet entering through the parameters Pq 
and Xg. In the limit Pq oc of infinite droplet viscosity, 
the function F(Xg,Xg) tends to zero, and we recover the 
result for fi obtained in Ref. [25]. 

In the weak dissipative limit of small boundary lay¬ 
ers, Sg,6g a, the dipole scattering coefficient for the 
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thermoviscous droplet reduces to 

3(1 + i) Po - 1 4 


_ 2 (po - 1 ) 
2 po + 1 


1 + 


1 + VQ’^f]Q ^ 2po + 1 a 
(Small-width boundary layers), 


(69a) 


bJ^ b:i"«1« bJ^ b:l^ b,|^ bill 


(69b) 


remarkable that for small particles suspended in a gas, 
where po ^ 1, the value of fi in Eq. (73) is three to 
five orders of magnitude larger than the value /i = 1 
predicted by ideal-fluid theory [3]. 


VII. RANGE OF VALIDITY 


2. A thermoelastic particle in a fluid 


In the long-wavelength limit Eq. (63) of a thermoelastic 
solid particle in a fluid, we obtain the result 


ft 


2 (p-Q - 1) (1 - G{x,)) 

2pQ -f 1 — 3G(xs) 


(70) 


with the function G(xs) given in Eq. ( 68 ). In this expres¬ 
sion, the only particle-related parameters are density and 
radius, and it is identical to that derived by Settnes and 
Bruns [25], who included the same two parameters in 
their study of scattering from a compressible particle in 
a viscous fluid using asymptotic matching. 

In the small-width boundary layer limit, (5g <C a, the 
dipole scattering coefficient for the thermoelastic solid 
particle /f reduces 


jsi _ 2 (po 


1 ) 


2 po + 1 


I + 3(1 + i) 


Po — I 4 


(7Ia) 




2pq + I a 
(Small-width boundary layers), 

(71b) 


which closely resembles Eq. (69) for /f. 


3. Asymptotic limits 


In the zero-width boundary layer limit, the dipole scat¬ 
tering coefficients /f and ff both reduce to the ideal- 
fluid expression [3], 


rfl ^ rsl _ 1) 

~ ~ 2po + 1 ’ 


(Zero-width boundary layers), 

(72) 


with the zero-width boundary layer limit defined for a 
droplet and a solid particle in Eqs. ( 6 Ib) and (67b), re¬ 
spectively. 

In the opposite limit of a point particle, = 

1/(2 - 1 - 3po) is finite and the expression for /q and /q is 
dominated by the G{xf) terms, with both cases yielding 
the asymptotic result 

fi = fi = ^(Po - 1), (Point-particle limit), (73) 

with the point-particle limit defined for a droplet and a 
solid particle in Eqs. (62b) and (65b), respectively. It is 


Before turning to experimentally relevant predictions 
derived from our theory, we discuss the range of validity 
of our results imposed by the three main assumptions: 
the time periodicity of the total acoustic fields, the per¬ 
turbation expansion of the acoustic fields, and the re¬ 
strictions associated with size, shape and motion of the 
suspended particle. 


A. Time periodicity 


The first fundamental assumption in our theory is the 
restriction to time-periodic total acoustic fields, which 
was used to obtain Eq. (3) for the acoustic radiation 
force evaluated at the static far-held surface dfli. Given 
a time-harmonic incident held, as studied in this work, 
a violation of time periodicity can only be caused by a 
non-zero time-averaged drift of the suspended particle. 
Denoting the speed of this drift by Up(t), we consider 
hrst the case of a steady particle drift. The assumption 
of time periodicity is then a good approximation if the 
displacement Ai is small compared to the particle radius 
a during one acoustic oscillation cycle r = 27r/uj used 
in the time averaging. A non-zero, acoustically-induced 
particle drift speed must be of second or higher order 
in 5ac, '^p/c ^ all hrst-order helds have a zero time 

average. Thus 


Ai v^T 27ru„ 27r Up A o 

- ^ -A- = -E =-P = _ 1 

a a uoa c a 


(74) 


and time periodicity is approximately uphe ld fo r reason¬ 
ably small perturbation strengths i/a/A, which is 

not a severe restriction in practice. In a given experi¬ 
mental situation, it is also easy to check if a measured 
non-zero drift velocity fulhlls UpT ^ a. 

In the case of an unsteady drift speed Up(t), the time- 
averaged rate of change of momentum in the huid 

volume bounded by dfti in Eq. (2) is non-zero, thus vio¬ 
lating the assumption (^) = 0 leading to Eq. (3). Only 
the unsteady growth of the viscous boundary layer in the 
fluid surrounding the accelerating particle contributes to 
(^), since equal amounts of momentum is fluxed into 
and out of the static fluid volume in the steady problem. 
Eor Eq. (3) to remain approximately valid, we must re¬ 
quire ( ^ ) to be much smaller than . To check this 
requirement, we consider a constant radiation force ac¬ 
celerating the particle. When including the added mass 
from the fluid, this leads to the well known time-scale 
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for the acceleration, 


2po + 1 


(75) 


Thus, small particles (a ^ S^) are accelerated to their 
steady velocity in a timescale much shorter than the 
acoustic oscillation period (rp ^ r), while the opposite 
(Tp ^ r) is the case for large particles (a ^ The un¬ 
steady momentum transfer to the fluid bounded by df^i 
is obtained from the unsteady part ^dmg(^) of fhe drag 
force on the particle as (^) = 7 / 0 ^'^dragW di- Using 
the explicit expression for Tdrag(^) given in Problem 7 
and 8 in §24 of Ref. [39], we obtain to leading order 


1 

-^rad 



2po + 1 n 


« 1 , 


2 a 

-- < 1 , 
TT 


for a ^ ^s, 
for a Sg. 


(76) 


We conclude that (^ ^ ^rad i^ both the large and the 
small particle limit, and hence the assumption of Eq. (3) 
is fulfilled in those limits. 

Considering typical microparticle acoustophoresis ex¬ 
periments, the unsteady acceleration takes place on a 
timescale between micro- and milli-seconds, much shorter 
than the time of a full trajectory. Typically, the unsteady 
part of the trajectory is not resolved and it is not impor¬ 
tant to the experimentally observed quasi-steady particle 
trajectory. In acoustic levitation [8-11], where there is 
no drift, the assumption of time periodicity is exact. We 
conclude that the assumption of time periodicity is not 
restricting practical applications of our theory. 


B. Perturbation expansion and linearity 

The second fundamental assumption of our theory is 
the validity of the perturbation expansion, which requires 
the acoustic perturbation parameter ^ac of Eq. (10) to 
be much smaller than unity. Eor applications in particle¬ 
handling in acoustophoretic microchips [13, 15], this con¬ 
straint is not very restrictive as typical resonant acoustic 
energy densities of 100 J/m^ result in ^ac ^ 10~^. 

Given the validity of the linear first-order equations, 
the solutions we have obtained for /q and fi based on the 
particular incident plane wave = (pQ are general, 
since any incident wave at frequency 00 can be written as 
a superposition of plane waves. 


C. Oscillations of the suspended particle 

The third fundamental assumption of our theory is the 
assumption of small particle oscillation amplitudes, al¬ 
lowing the boundary conditions to be evaluated at the 
fixed interface position r = a. This assumption puts 


physical constraints on the volume oscillations, Eig. 1(a) 
and (b), and the center-of-mass oscillations, Eig. 1(c). 

The volume oscillations of the particle are due to me¬ 
chanical and thermal expansion. Erom the definition of 
the compressibility Kg and the volumetric thermal expan¬ 
sion coefficient we estimate the maximum relative 
change in particle radius Aa/a to be 


Aa 

a 

Aa 

a 


k' 

- -fpi = y^ac < 1, 

- yUi - ^ (7 - l)ap£ac < 1- 


(77a) 

(77b) 


Here, we have used and ob- 

tained from Eq. (14) in the adiabatic limit = 0 com¬ 
bined with Eq. (16). Except for gas bubbles in liquids, 
for which ^ 1 , these inequalities are always fulfilled 
for small perturbation parameters ^^c- 

The velocity of the center-of-mass oscillations is found 
from Eq. (37) of Ref. [25] to be = f f^e 

large-particle limit, fi is given by Eq. (72), which implies 
0 < < 3r’in: where the lower and the upper limit is 

for po ^ 1 and p ^ 1, respectively. In the point-particle 
limit, Eq. (73), independent of pg. The relative 

displacement amplitude Ai/a is hence estimated as 


A /) osc 

A£ ^ Up 
a Lua 


3 A 
2po + 1 21 : 

A 


« 1 , 

« 1 , 


for a ^ Sg, 
for a Sg, 


(78) 


and thus the general requirement is that ^ ^'KajX. 
Eor large particles in typical experiments, this restriction 
is not severe. However, for small particles it can be re¬ 
strictive. Eor example, to obtain Aija < 0.05, we find 
for particles of radius a = 100 nm in water at 1 MHz 
and particles of radius a = 1 pm in air at 1 kHz, that 
^ac ^ 10 ”^ and 5ac < 10 “^, respectively. 


VIII. MICROPARTICLES AND DROPLETS IN 
STANDING PLANE WAVES 


The special case of a one-dimensional (ID) standing 
plane wave is widely used in practical applications of the 
acoustic radiation force in microchannel resonators [ 12 - 
24] and acoustic levitators [8-11]. The many application 
examples as well as its relative simplicity, makes the ID 
case an obvious and useful testing ground of our the¬ 
ory. In the following, we illustrate the main differences 
between our full thermoviscous treatment and the more 
conventional ideal-fluid or viscous-fluid models using the 
typical parameter values listed in Table H. 

We consider a standing plane wave of the form = 
PaCOs(fcoy), Vin = ^ Pa sin(fco 2 /)ey, with acoustic energy 
density = jK^pl = jPqvI, where p^ and are the 
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TABLE IL Material parameter values at ambient pressure 0.1 MPa and temperature 300 K used in this study, given for water 
(wa) [40-43], an average liquid food oil [44], air [45], and polystyrene (ps) [46-49]. Parameter values for water and oil at other 
temperatures are obtained from the fits in Refs. [40, 44]. 

Parameter Symbol Value (wa) Value (oil) Value (air) Value (ps) Unit 


Longitudinal speed of sound 

c 

1.502 

X 

10^ 

1.445 

X 

10° 

3.474 

X 

10^ 

2.40 

X 

10° 

-1 

m s 

Transverse speed of sound 

Ct 


- 



- 



- 


1.15 

X 

10° 

-1 

m s 

Mass density 

Po 

9.966 

X 

10^ 

9.226 

X 

10^ 

1.161 

X 

10° 

1.05 

X 

10° 

kg m“° 

Compressibility 

Kg 

4.451 

X 

10“^° 

5.192 

X 

10“^° 

7.137 

X 

10“° 

2.38 

X 

10“^° 

Pa“^ 

Thermal expansion coefficient 

OLp 

2.748 

X 

10“^ 

7.046 

X 

10“^ 

3.345 

X 

10“° 

2.09 

X 

10“^ 

K“^ 

Specific heat capacity 

Cp 

4.181 

X 

10^ 

2.058 

X 

10° 

1.007 

X 

10° 

1.22 

X 

10° 

Jkg“^ K“^ 

Heat capacity ratio 

7 

1.012 

X 

10° 

1.151 

X 

10° 

1.402 

X 

10° 

1.04 

X 

10° 

1 

Shear viscosity 

Vo 

8.538 

X 

10-^ 

4.153 

X 

10“^ 

1.854 

X 

10“° 




Pa s 

Bulk viscosity ^ 

b 

Vo 

2.4 

X 

10“° 

8.3 

X 

10“^ 

1.1 

X 

10“° 




Pa s 

Thermal conductivity 

Pth 

6.095 

X 

10“^ 

1.660 

X 

10“^ 

2.638 

X 

10“^ 

1.54 

X 

10“^ 

Wm“^ K“^ 


^ The bulk viscosity is negligible for scattering in the long-wavelength limit but has been included for completeness. Values for water, oil 
and air are estimated from Refs. [50], [51], and [52], respectively. For oil, 77 q is obtained from the attenuation constant olq at 298.15 K 
and 10 MHz [51] using oq = 27r^/^/(PoC^)[^o + + (t - l)^th/cp]. 


pressure and the velocity amplitude, respectively. Ex¬ 
pression (5) for the radiation force then simplifies to 

r’i'D* = 47r i'aca^fco-E’ac sin(2fcoy)ej,, (79a) 
^ac = [/o] + f Re [/i], (79b) 

where is the so-called acoustic contrast factor. The 
radiation force is thus proportional to which con¬ 
tains the effects of thermoviscous scattering in /q and fi. 
Note that for positive acoustic contrast factors, > 0? 
the force is directed towards the pressure nodes of the 
standing wave, while for negative acoustic contrast fac¬ 
tors, < 0, it is directed towards the anti-nodes. 

The acoustic contrast factor evaluated di¬ 

rectly for an arbitrary particle size by using the expres¬ 
sions for the scattering coefficients, either /q and /f for 
a fluid droplet or /q and fi for a solid particle. For 
ease of comparison to the work of King [1], Yosioka and 
Kawasima [2], and Doinikov [4-6], we give the expression 
for the acoustic contrast factor <4>ac of a fluid droplet for 
small boundary layers and in the point-particle limit. In 
the small-width boundary layer limit one obtains 


incompressible particles for which = 0. The second 
term is the viscous correction, which agrees with the re¬ 
sult of Settnes and Bruns [25] for infinite particle vis¬ 
cosities, but extends it to finite particle viscosities. Note 
that the viscous correction yields a positive contribution 
to the acoustic contrast factor, while the thermal correc¬ 
tion from the third term is negative. The result given in 
Eq. (80) is in agreement with the expression for the radia¬ 
tion force in a standing plane wave given by Doinikov [6] 
in the weak dissipative limit of small boundary layers. 
However, this is not seen without considerable effort com¬ 
bining and reducing a number of equations. Although we 
find Doinikov’s approach rigorous, it lacks transparency 
and is difficult to apply with confidence. 

In the point-particle limit of infinitely large boundary 
layer thicknesses compared to the particle size, we obtain 


Kc = 1 (1-kJ - (l-Po) - (7-1 )PoCp(i-tY) 

3 [ poCp/ 

(81a) 

(Point-particle limit), 

kcit Wc? 1 and lx'1^ \x^\^, |x(|^ 0, (81b) 


l/^_ 

3\2po + l 


+ 


l + Vo 


Po - 1 

2 po + 1 


7-1 


21 + 


Ah Ah 


1 - 


Po^p J 


7 A 


(80a) 


(Small-width boundary layers), 
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(80b) 


The first term is the well-known result given by Yosioka 
and Kawasima [2], which reduces to that of King [1] for 


in agreement with the viscous result of Settnes and 
Bruns [25], when omitting the last term stemming from 
thermal effects. The result for Eq. (81) is written 

in a form which emphasizes how parameter contrasts be¬ 
tween particle and fluid lead to scattering. As expected, 
for kg = 1 and pQ = the scattering due to compress¬ 
ibility and density (inertia) mechanisms vanishes. This is 
true for large particles [1-3, 25], and it is reasonable that 
it remains true in the point-particle limit. The expres¬ 
sions for the acoustic radiation force on a point-particle 
in a standing plane wave given by Doinikov [4-6] do not 
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FIG. 3. (Color online) Acoustic contrast factor plotted as a function of (5g/a, the viscous boundary layer thickness in 
the medium normalized by particle radius. The curves are calculated using ideal theory (green), viscous theory (blue), and 
thermoviscous theory (red) for (a) an oil droplet in water (wa-oil) and (b) a water droplet in oil (oil-wa), both at 20° C. 
Thermoviscous theory leads to corrections to the acoustic radiation force around 100%. The vertical dashed lines indicate 
examples of particle sizes corresponding to the given value of S^/a at / = IMHz. Note that the acoustic contrast factor changes 
sign at a critical particle radius for the case of water droplets in oil considered in (b). 


have this property, which is likely due to a sign-error or 
a misprint in the term corresponding to our dipole scat¬ 
tering coefficient fi in the point-particle limit Eq. (73), 
as was also suggested by Settnes and Bruus [25]. 

The small-width boundary layer limit and the point- 
particle limit are useful for analyzing consequences of 
thermoviscous scattering on the acoustic radiation force, 
but we emphasize that our theory is not restricted to 
these limits. In general, the scattering coefficients /q and 
fi are functions of the non-dimensionalized wavenumbers 
Xg, Xt, Xg, and x[ . These may all be expressed in terms of 
the particle radius a normalized by the thickness of the 
viscous boundary layer in the medium (5g, 


droplet (oil-wa) using the parameters of a typical food 
oil given in Table II. Since the density contrast of wa¬ 
ter and oil is small, the dipole scattering with its viscous 
effects is small, while on the other hand the thermal ef¬ 
fects in the monopole scattering are significant. This is 
clearly seen from Fig. 3, where the acoustic contrast fac¬ 
tor is plotted for the two cases as function of d^/a 
using ideal theory, viscous theory, and full thermoviscous 
theory. Fig. 3 shows that for sub-micrometer droplets at 
MHz frequency the thermoviscous theory leads to cor¬ 
rections around 100% as compared to the ideal and the 
viscous theory, which manifestly demonstrates the im¬ 
portance of thermal effects in such systems. 


_ /I , ^ ' 

Xg (1 T i) , Xg 



(82a) 

(82b) 


where we have used i5s = Sg^^fj^/pQ, (5t = S^yX/Pv, 
K = 4^[(1 -^')Ah]/Pr, with Pr = J/o/Ah being the 

Prandtl number of the fluid medium and set to zero 
for the fluid droplet case. Below, we investigate the ther¬ 
moviscous effects on the acoustic radiation force by plot¬ 
ting the acoustic contrast factor ^ function of S^/a, 

ranging from zero boundary-layer effects at (5g/a = 0 to 
maximum effects in the limit S^/a oo. 


A. Oil droplets in water and water droplets in oil 

We first consider the cases of water with a suspended 
oil droplet (wa-oil) and of oil with a suspended water 


We note from Fig. 3(a) that the acoustic contrast fac¬ 
tor of oil droplets in water is negative, which means that 
oil droplets are focused at the pressure anti-nodes. Con¬ 
versely, water droplets in oil are thus expected to be fo¬ 
cused at the pressure nodes. However, in Fig. 3(b) we 
see that thermoviscous theory predicts a tunable sign- 
change in the acoustic contrast factor as a result of the 
negative thermal corrections to the monopole scattering 
coefficient. This means that droplets above a critical size 
threshold experience a force directed towards the pres¬ 
sure nodes, while droplets smaller than the threshold ex¬ 
perience a force towards the anti-nodes, even though the 
only distinction between the droplets is their size. This 
sign-change in can also be achieved for elastic solid 
particles under properly tuned conditions. By changing, 
for example, the compressibility contrast the curves 
for ^ac(^s/^) be shifted vertically and a possible size- 
threshol d conditio n may be changed. Moreover, since 
4 = VAo/(Po‘^) and <5^ = y/2k^h/(poCpU;), there are 
several direct ways of tuning a threshold value, e.g. by 
frequency or by changing the density of the medium. 
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FIG. 4. (Color online) Acoustic contrast factor 4)^^, for particles in air plotted as a function of 8^/a^ the viscous boundary layer 
thickness in the medium normalized by particle radius. The curves are calculated using ideal theory (green), viscous theory 
(blue), and thermoviscous theory (red) for (a) a polystyrene particle in air (air-ps) and (b) a water droplet in air (air-wa), both 
at 300 K. Ideal theory predicts a constant value of 4>ac = 5/6 independent of particle size. For particles much smaller than the 
boundary layer thickness, however, thermoviscous theory predicts huge deviations from ideal theory leading to acoustic contrast 
factors two orders of magnitude larger than expected from ideal-fluid theory. The vertical dashed lines indicate examples of 
particle sizes corresponding to the given value of (5s/a at /= 1 kHz. 


B. Polystyrene particles and water droplets in air 

Using the particular cases of a polystyrene particle and 
a water droplet suspended in air as main examples, we 
study the effects of a large density contrast po ^ 1, for 
which our thermoviscous theory predicts much larger ra¬ 
diation forces on small particles than ideal-fluid theory, 
for which = 5/6 independent of particle size. This 

is demonstrated in Fig. 4, where ^ac is plotted as a func¬ 
tion of 5^/a for the two particle types. In the large- 
particle limit (5s/a = 0 , boundary-layer effects are neg¬ 
ligible, and ideal, viscous, and thermoviscous theory pre¬ 
dict the same contrast factor <Fac = 5/6, but as d^/a in¬ 
creases, the thermoviscous and viscous theory predict an 
increased value of <Fac, approximately 2 <Fac^^^ for ^g/a = 1 
as seen in the insets of Fig. 4(a) and (b). Decreasing the 
particle size further, 5^/a ^ 1 , the thermoviscous effects 
become more pronounced with <Fac/^Lc ^ ^ 10^- Choos¬ 
ing the frequency to be 1 kHz, this remarkable deviation 
from ideal-fluid theory is obtained for moderately-sized 
particles of radius a ~ 2 pm. 

While in Fig. 4(a) for the polystyrene parti¬ 

cle is a monotonically increasing function of (5s/a, the 
in Fig. 4(b) of a water droplet exhibits a non¬ 
monotonic behavior. For small values of 5^/a < 25, the 
viscous dipole scattering dominates resulting in a pos¬ 
itive contrast factor < 10^. For larger values, 

5^/a > 25, thermal effects in the monopole scattering be¬ 
come dominant leading to a sign-change in and 

finally to lar|e negative contrast factors approximately 
equal to — 10 "^ as the point-particle limit 5^/a ^ 1 is ap¬ 
proached. This example clearly demonstrates how the 
acoustic contrast factor may exhibit a non-trivial size- 


dependency with profound consequences for the acoustic 
radiation force on small particles. The detailed behavior 
depends on the specific materials but can be calculated 
using Eq. (79) and the expressions for /q and fi listed in 
Table I. 


IX. CONCLUSION 

Since the nominal work of Epstein and Carhart [26] 
and Allegra and Hawley [27], effects of thermoviscous 
scattering have been known to be important for ultra¬ 
sound attenuation in emulsions and suspensions of small 
particles. In this paper, we have by theoretical analysis 
shown that thermoviscous effects are equally important 
for the acoustic radiation force on a small particle. 

is evaluated from Eq. (5), or more generally from 
Eq. ( 6 ), using our new analytical results for the ther¬ 
moviscous scattering coefficients /q and fi summarized 
in Table I. Our analysis places no restrictions on the 
viscous and thermal boundary layer thicknesses 6^ and 
(5t relative to the particle radius a, a point which is es¬ 
sential to calculation of the acoustic radiation force on 
micro- and nanometer-sized particles. 

The discussion in Section II leading to Eq. (5) for 
as well as the discussion of the range of validity presented 
in Section VII, are intended to provide clarification and 
a deeper insight into the fundamental assumptions of the 
theory for the acoustic radiation force. Eoremost, we 
have extended the discussions of the role of streaming, 
the fundamental assumption of time periodicity, and the 
trick of evaluating the radiation force in the far-held. To 
our knowledge, the exact non-perturbative expression (3) 
for the radiation force evaluated in the far-held has not 
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previously been given in the literature. 

For the simple case of a ID standing plane wave at a 
single frequency, the expression (6) for simplifies to 
the useful expression given in Eq. (79), which involves 
the acoustic contrast factor Similar simplified ex¬ 

pressions can be derived for other cases of interest such 
as that of a ID traveling plane wave. An important result 
from the discussion of the simple ID case in Section VIII 
is that we must abandon the notion of a purely material- 
dependent acoustic contrast factor general, <Fac 

also depends on the particle size, and in many cases this 
size-dependency can even lead to a sign change in 
a critical threshold. Recent acoustophoretic experiments 
on sub-micrometer-sized water droplets and smoke par¬ 
ticles in air may provide the first evidence of this predic¬ 
tion [53]. Considering only viscous corrections, however, 
the authors could not fully explain their data. Our analy¬ 
sis suggests that thermoviscous effects must be taken into 
account when designing and analyzing such experiments. 

Our results for the acoustic radiation force in a stand¬ 
ing plane wave evaluated using Eq. (79) agree with the 
expressions obtained from the work of Doinikov [4-6] in 
the limit of small boundary layers, but not in the opposite 
limit of a point particle. In our theory both of these lim¬ 
its are evaluated directly using the derived analytical ex¬ 
pressions valid for arbitrary boundary layer thicknesses, 
and we have furthermore given a physical argument sup¬ 
porting our result in the point-particle limit. Considering 
the viscous theory of Danilov and Mironov [7], we remark 
that their result is based on the viscous reaction force on 
an oscillating rigid sphere [39] instead of a direct solution 
of the governing equations for an acoustic field scattering 
on a sphere. 

Importantly, we have shown that the acoustic radiation 
force on a small particle including thermoviscous effects 
may deviate by orders of magnitude from the predictions 
of ideal-fluid theory when there is a large density contrast 
between the particle and the fluid. This result is particu¬ 
larly relevant for acoustic levitation and manipulation of 
small particles in gases [8-11]. Thermoviscous effects can 
also be significant in many lab-on-a-chip applications in¬ 
volving ultrasound handling of submicrometer-sized par¬ 
ticles such as bacteria and vira [18, 20]. 

A firm theoretical understanding of thermoviscous ef¬ 
fects, and of the particle-size-dependent sign change of 
the acoustic contrast factor, could prove important for 
future applications relying on ultrasound manipulation 
of micro- and nanometer-sized particles. 


Appendix A: Velocity and normal stress in spherical 
coordinates 


In spherical coordinates (r, 6>, ip) with azimuthal sym¬ 
metry, using that i;i = V(/) + Vx'0 with (p = (p^ cj)^ 


and ip = pjg the first-order velocity components are 

Vir = dr4>^ -(Ala) 

Vw = \dg(i> - [r V’s] • (Alb) 

Inserting this into Eq. (52), we obtain the normal com¬ 
ponents of the first-order stress tensor 


C^lrr = %(2fcc - (f) 


sin^ 


sin l9 ( - 4^ 

r T 


, (A2a) 


1 


o-ier = ‘^Vodg V’s - 


+ —Oq 

r 


I 


sinO 


a0(sin6» V>s) 


(A2b) 


Appendix B: The scattering coefficients /o and /i 


Here, we outline the calculation of the monopole and 
dipole scattering coefficients /q and fi in the long- 
wavelength limit where the particle radius and the 
boundary layer thicknesses are assumed much smaller 
than the wavelength. Defining the small parameter 
5 = kQa I, we note that 

a fluid particle furthermore are all of order 5. The 
calculation is carried out to first order in 5. 


1. The monopole scattering coefficient /q 


The monopole scattering coefficient /q may be ob¬ 
tained from Eqs. (56a), (56c), (56d) and (56f) setting 
n = 0 and Cq = Cq = 0. All Bessel functions of the 
small arguments ^ s ^ 1 are expanded to first 

order in 5 using Eq. (C5) of Appendix C, and in the (un¬ 
primed) fluid medium we neglect xl in comparison to x ‘^. 
Thus, we arrive at 

Aq— + A!q\x'c - BQXth^{x^) + BQx[ji{x[) = be, 

X Q O O 

(Bla) 

Ao^c (^1 - - Ao6e + Bobtho{xt) - Bob'tjo(x[) = -b^, 

(Bib) 

AhhK— + Ao-fc(h^cb - -Bo^th&ta^t^i(a:t) 

Xc 3 

+ Bokthb[x'jj^{x[) = ^Ki^b^xl, (Blc) 
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^oVo 


(4-4)— +a;s 


^oVo 


/2 


+ Borio {xi-2xt)ho{x^)-2xlho{x^) 


-BWo 


(x' - 2x[ )io(x[) - 2xfio (x( 


= -VqXs , 

(Bid) 


where Eq. (C3) is used to write ^o(^) = ~9i{^) 
spherical Bessel of Hankel function go{x). 

Multiplying Eq. (Blc) by 1/(/^th^t) using the ratios 


b X 

J^ = -(7-l)% 

°t Xt 




1 


bi 


bt Xdp-Oth 


Uq L,p 

bt bt bf. 


(B2) 




Cp 


2 

c 

2 ’ 


of the 6 -coefficients defined in Eq. (50) (here, y = 1 for 
a droplet and X = X ^ solid particle, respectively, 
while Eqs. (16), ( 22 ), and (39e) is used to reduce h'^jh^)^ 
we note that the Aq and A'q terms can be neglected to 
order 5, and we obtain 

_ /Cth&t x^h^{x^) ^ , dp xXi{xt) ^ 

^0 — TTTT / . / /N ^0 — X T-r ~T~T1T ^ 0 - 


Khb't 4ii(4) 


PoCp x(ji'i(4) 


With this, we eliminate Bq from the system of equa¬ 
tions (Bl), and the remaining three equations become 


t 


1 '2 

3^c 

-5i\ 

/Ao\ 


/ 2 . 

/£c \ 

3 

U (j_ . 

6t Vaic 

-1) 

K 

bt 

S2 

-2 

X 

= 

be 

bt 
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where we have introduced the functions Si, S 2 , and S^, 


iS-functions are expressed in terms of the dimensionless 
wavenumbers as 
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where H{xi,x[) is given in Eq. (59b). The coefficient 
Aq is now found from Eq. (B4) by the method of de¬ 
terminants (Cramer’s rule) as Aq = D{Aq)/D, where 
D is the determinant of the left-hand-side system ma¬ 
trix and D{Aq) is determinant of the system matrix in 
which the first column (the Aq coefficients) are replaced 
by the right-hand-side column with the inhomogeneous 
terms. The monopole scattering coefficient /q in the long- 
wavelength limit can then be expressed as 


. _ 3i , _ 3i D{Ao) 

Jo — 3 ^ 0 — 3 y-v 5 

with the determinants D and D{Aq) given by 
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(B5c) 


and the relative shear constant 7)0 obtained from 
Eq. (41b), 


(B9b) 


The solution Aq = D{Aq)/D, though written somewhat 
differently, agrees with Allegra and Hawley’s Eq. (10) of 
Ref. [27]. 


~ _1lQ _ ~ X^ 
Vo - — - Po —. 
Vo X, 


(B 6 ) 


In obtaining the expression for Sq we have used Eq. (C3) 
to substitute gQ{x) = —gQ{x)-\-{2/x)gi{x) for any spheri¬ 
cal Bessel or Hankel function g{x). Using Eq. (B 2 ), (B 6 ), 
and the explicit forms (C4) of the Bessel functions, the 


/o fo'B a suspended thermoviseous droplet 

Eor a suspended thermoviseous droplet, the precise def¬ 
inition of the long-wavelength limit is given in Eq. (58). 
In this case, the shear mode characterized by x^ inside 
the droplet corresponds to a boundary layer, and con¬ 
sequently comparison to the compressional mode inside 
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and outside the droplet yields x^jx'^ ^ I ^ ^1. 
This, combined with bc/b^ r\j b'c/bt ~ xl/xl ~ < 1 
from Eq. (B2), leads to the following simplification of 
Eq. (B9) to first order in e, 


D - - 


D{Ao) - 


(BlOa) 



PoCp ^ 

(BlOb) 


When inserting this into Eq. (B8), we obtain 

/o« = !-«, +3(7(Bll) 

V PoS / *82 

which upon substitution with ^ = (l — -^-^H{x^^x[) 
from Eq. (B7) with y = 1, leads to the final analytical 
result for /q given in Eq. (59). 


b- fo for CL suspended thermoelastie partiele 


layer mode to a propagating transverse wave with x'^ ^ 
s ‘^. A further implication is that the shear constant ratio 
of Eq. (B6) becomes large, fjQ = pox^/x'^ ^ ^ 1, 

and order of magnitude wise, the ^'-functions of Eq. (B7) 
obey Si ^ S 2 ^ Combining this with the following 

expression derived from Eqs. (39e), (48), and (B6), 

Xc -^s I = -XPo^s, (B12) 



the leading-order expansions in 5 of the determinants D 
and D{Aq) in Eq. (B9) become 
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The qualitative change going from the thermoviscous 
droplet to the thermoelastie particle lies in the shear 
mode, which changes from a highly damped boundary 


Erom this and Eq. (B8), we obtain the monopole scatter¬ 
ing coefficient /q for a thermoelastie particle suspended 
in a thermoviscous fluid. 
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Erom Eq. (B7) we obtain the leading-order expansions in 
5 for the ratios Si/S 2 and S^/S 2 , 
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(B15) 


with the function H{xi,x'i) defined in Eq. (59b). Insert¬ 
ing this into Eq. (B14) and using Eqs. (B2) and (B6), 
and the expression (39e) for x 1 we arrive at the final 
analytical form for /q given in Eq. (64). 


2. The dipole scattering coefficient fi 

In the long-wavelength limit, for each order n, the 
terms containing and and thus the variables x^ 
and x[, in the system of boundary equations (56) are of 
negligible order relative to the terms containing 


and the inhomogeneous terms. Eormally, this is 
seen by writing up and inverting the entire 6-by-6 matrix 
equation for the six coefficients for a given n. A quicker 
way to see this, is to write Eqs. (56c) and (56d) as 


f ^n(^t) A f A 

V ) \ K ) 

f ^n^n(^c) Jni^c) A 

\ ^n^n(^c) Jni^c) J 


(B16) 


where we have used ^ ^ and ^^ 1. Insert- 

h ^ h bt^ /cth 

ing the expressions for B^ and Br^ obtained by inversion 
of this equation into Eqs. (56a), (56b), (56e), and (56f), 
we see that due to the factor s‘^ each term related to B^ 
or B'^ are negligible in all four equations. In treating 
Eq. (56e) it might be useful to use the Bessel’s equa¬ 
tion (C2). Consequently, returning to the dipole prob¬ 
lem with n = 1, terms with Bi^B[ are omitted and the 
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system of equations reduces to four equations with four 
unknowns, namely Eq. (56a), Eq. (56b), Eq. (56e), and 
Eq. (56f) without the terms of Eor n = 1 we thus 

obtain the simplified system of equations 


Xcfiixc) + A:^x^h[{x^) - 

ji{x^) + - Cl [a;g/ii(a;s) + /ii(a;s)] 

= Aj^{x^) - C[ [xJi{x^) + ji(a;')] , 


Vo 


Vo 


XckiXc) + Aix^h2(x^) + ^Cixlh'lix^) 


= Vo 


A'lXcj2{x'c) + -^C[xl^j”{x'^) 


XsJliXc) - ^xj2{Xc)\ - 4:C-i^VoXsh2{Xs 
+ Aiijo xlhi{xc) - Ax^h2{x^) 

= AWo x'^jiix'^) - 4xj2iXc) 

- 4:C[riox'j2{Xs), 


(B17a) 

(B17b) 

(B17c) 


(B17d) 


where we have rewritten the last two equations using the 
recurrence relations obtained from Eq. (C3) 


relevant scalings in the long-wavelength limit for the fluid 
droplet and the solid particle, respectively. 

fi for a suspended thermoviseous droplet 

In the long-wavelength limit for the fluid droplet case 
the scalings of Eq. (58) apply. Using the approximate 
expressions for the spherical Bessel and Hankel functions 
Eq. (C5) applicable for small arguments and examining 
the resulting system of equations (B19) one finds that 
some terms may be omitted to first order in 5 . The sim¬ 
plified system of equations (B19) for the fluid droplet 
case takes the form 

--^Ai + C-yX^h 2 {x^) - C[xj 2 {Xs) = 0, (B20a) 

Xc 


-2CiX^ho(x^) - + 2C[x'jo{x^) = -x^, (B20b) 

-^Ai + ^Cix^h'lix^) - ^C[fioxfji{Xs) = 0, (B20c) 

Xr. ^ ^ 


xg'i (x) - 5 i (x) = -xg 2 (x), (B18a) 

g'lix) =-giix)+ ^g 2 ix), (B18b) 


—+ 6Ci/ii(xJ + A[pox'^ - QC[poji{x's) = x^, 

Xc 

(B20d) 


valid for any spherical Bessel or Hankel function g. 

Simplifying the system of equations we multiply 
Eq. (B17a) by (—1) and add to it Eq. (B17b), then use the 
recurrence relation (B18a). Eq. (B17b) is multiplied by 2 
and Eq. (B17a) is added while using the recurrence rela¬ 
tion xg'i{x) ^2gi{x) = xgQ{x). We leave Eq. (B17c) as it 
is. To Eq. (B17d) we add 4 times Eq. (B17c) and use the 
recurrence relation (B18b). With some rearrangements, 
these manipulations give 


^ 1 X 0 / 12 ( 2 ^ 0 ) + CiXs/l2(Xg) 

- A[x'j 2 ixc) - C[x'j 2 {x's) = -xj 2 {xc): (B19a) 
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(B19c) 


AiK{x^) - 2 Ci/ii(xg) 

- Po [A'ljiix'c) - 2C(ji(a;s)] = -3i{Xc)- (B19d) 


where fj^x^ = p^xl was used to simplify the last equa¬ 
tion. The equations may be further simplified using the 


Subtracting Eq. (B20c) from Eq. (B20a) and using 
Eq. (B18b), we can express C[ by Ci, 


^ VoXsQiXs) 

Q{Xs) = x'sjlix'^) - 2(1 - ^)j2{Xs)- 


(B21a) 

(B21b) 


Then, using this relation to eliminate C[ in Eq. (B20a), 
we arrive at the first of the two equations in Eq. (B22). 
The second equation (B22b) is obtained by adding 
Eq. (B20b) and Eq. (B20d) in order to eliminate then 
making use of the recurrence relation 3gi{x) — xgQ{x) = 
xg 2 {x). The resulting two equations for Ai and Ci are 


Xc 
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From this, and using again the relation 3^i(x) — xgQ{x) = xg 2 {x), we obtain the dipole expansion coefficient 




_ ^xljpo - l)[h2{Xs)rjoQ{x'^) - x^hi{x^)j2{x'^)] _ 

[3/i2(a;s) - 2(po - l)^o(a;s)]%Q(4) - (2po + l)a;s/ii(a;s)j2(4)' 


(B23) 


This result, but with a small error in the numerator, was 
first obtained by Epstein and Carhart [26]. We reduce 
the fraction by ^oQ(^s)^o(^s) explicit ex¬ 

pressions for the Bessel and Hankel functions in Eq. (C4) 
to introduce the functions G{Xg) and F{xg,x'g) given ex¬ 
plicitly in Eqs. (68b) and (68c), respectively. 


Appendix C: Special functions 

The Legendre differential equation solved by Legendre 
polynomials (cos 0) of order n is [54] 


G(xs) — 1 + 


h2iXs) 

^0 (^s) 


flohoiXs)Qix:)' 


(B24a) 

(B24b) 


Then, using that fi = —hia;^. we arrive at the final 
expression (68a) for the dipole scattering coefficient /f. 


si^^ (^sin6»Tp^(cos0)^ + n(n + l)P„(cos0) = 0. 

(Cl) 

The Bessel differential equation solved by spherical 
Bessel or Hankel functions gn{x) of order n is [54] 

x"^ [g'nix) + gn(x)] = n(n + l)gn(x) - 2xgn(x), (C2) 


fi for a suspended thermoelastie partiele 


with a prime indicating differentiation with respect to 
the argument. Useful recurrence relations for gn{x) are 


In the long-wavelength limit for the solid particle the 
scalings of Eq. (63) apply. Using the approximate ex¬ 
pressions for the spherical Bessel and Hankel functions 
Eq. (C5) applicable for small arguments and examining 
the resulting system of equations (B19) one finds that 
some terms may be omitted to first order in 5 . The sim¬ 
plified system of equations (B19) in the solid particle case 
takes the form 


A ^ 

dx 

A ^ 

dx 


n+l 


9n (^) 
gn{x) 


= -X "5n+i(a;), 
= X^^'^Qn-iix). 


(C3a) 

(C3b) 


The lowest-order spherical Bessel functions jn{x) and 
Hankel functions of the first kind h^{x) are [54] 


- 2 + C'i^s^ 2 (^s) = 0: (B25a) 
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(B25c) 

- - 6Gihi(xJ - A[pqx'^ + 2Gipo^s = 

(B25d) 


Multiplying Eq. (B25b) by (—Po) and adding it to 
Eq. (B25d), then substituting Ci using Eq. (B25a), and 
finally using the recurrence relation 3pi(x) — xgQ(x) = 
xg 2 {x)^ leads to the expansion coefficient A^, 


Ai — 


3 ^c(Po ^)^ 2 (^s) 

3h2(xs) - 2{po - l)ho{xs)' 


(B26) 


Again, using that fi = —Qix^ ^ Ai and introducing 
G(Xg) as defined in Eq. (B24a), we obtain after some 
rearrangement the final result for given in Eq. (70). 
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